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ABSTRACT 



We present the results of a systematic study of the formation and evolution of bi- 
naries containing black holes and normal-star companions with a wide range of masses. 
We first reexamine the standard formation scenario for close black-hole binaries, where 
the progenitor system, a binary with at least one massive component, experienced a 
common-envelope phase and where the spiral-in of the companion in the envelope 
of the massive star caused the ejection of the envelope. We estimate the formation 
rates for different companion masses and different assumptions about the common- 
envelope structure and other model parameters. We find that black-hole binaries with 
intermediate- and high-mass secondaries can form for a wide range of assumptions, 
while black-hole binaries with low-mass secondaries can only form with apparently 
unrealistic assumptions (in agreement with previous studies). 

We then present detailed binary evolution sequences for black-hole binaries with 
secondaries of 2 to 17 Mq and demonstrate that in these systems the black hole can 
accrete appreciably even if accretion is Eddington limited (up to 7 Mq for an initial 
black- hole mass of 10 M©) and that the black holes can be spun up significantly in the 
process. We discuss the implications of these calculations for well-studied black-hole 
binaries (in particular GRS 1915-1-105) and ultra-luminous X-ray sources of which GRS 
1915-1-105 appears to represent a typical Galactic counterpart. We also present a de- 
tailed evolutionary model for Cygnus X-1, a massive black-hole binary, which suggests 
that at present the system is most likely in a wind mass-transfer phase following an 
earlier Roche-lobe overflow phase. Finally, we discuss how some of the assumptions in 
the standard model could be relaxed to allow the formation of low-mass, short-period 
black-hole binaries which appear to be very abundant in Nature. 

Key words: binaries: close - stars: black holes - X-rays: stars - stars: individual: 
GRS 1915-1-105 - stars: individual: Cygnus X-1 - gravitation 



1 INTRODUCTION 

There are currently 17 binary systems containing black holes 
for which dynamical mass estimates are available (see e.g. 
Table 1 of Lee, Brown & Wijers 2002 [LBW], and references 
therein; Orosz et al. 2002). According to conventional wis- 
dom, these systems formed from primordial binaries where 
at least one of the stars was quite massive (i.e. M>20- 
25 M0). If mass transfer from the primary to the secondary 
commences at an orbital period in the range of ~ 1-10 yr, 
a common envelope may form during which the hydrogen- 
rich envelope of the primary is expelled (Paczyhski 1976) . If 
the secondary and the core of the primary avoid a merger, 
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then the massive core may evolve to core collapse and the 
formation of a black hole in a close binary. 

For 9 of the 17 black-hole binaries (see e.g. LBW), the 
current-epoch companion mass is < 1 Mq and the orbital 
periods are < 1 day. For reasons discussed later in the text, 
these systems probably had primordial secondaries whose 
mass was not substantially greater than ~ 1.5 Mq (but see 
§ 4.4). One quantitative difficulty with the common-envelope 
scenario for forming this type of black-hole binary is that the 
amount of orbital energy that can be released by the spiral- 
in of a low-mass secondary may not be sufficient to eject the 
massive envelope of the primary. It has long be recognized 
that this is energetically challenging even if the common- 
envelope ejection mechanism is very efficient (Podsiadlowski, 
Cannon & Rees 1995; Portegies Zwart, Verbunt & Ergma 
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1997; Kalogcra 1999; see, however, also Romani 1992). Fur- 
thermore, recent determinations of the binding energy of the 
envelopes of massive supergiants by Dewi & Tauris (2000, 
2001) suggest that all studies so far may have significantly 
underestimated how tightly bound these envelopes actually 
are, which seriously aggravates the problem. On the other 
hand, it has been estimated that there may be up to several 
thousand low-mass black-hole transients in the Galaxy (Wi- 
jers 1996; Romani 1998). This has led to several alternative 
formation scenarios for low-mass black-hole binaries, where 
either the low-mass companion is a third star in a triple 
system being captured into a tight orbit when the two mas- 
sive components merge (Eggleton & Verbunt 1986), or where 
the low-mass star forms after the black hole - out of a col- 
lapsed massive envelope (Podsiadlowski et al. 1995). In the 
present work we reexamine the standard formation scenarios 
for low-mass black hole binaries with plausible modifications 
to some of the usual assumptions 

By contrast, for 4 of the black-hole binaries, the mass 
of the companion is substantially larger (i.e. >6Mq), and 
the availability of orbital binding energy for ejecting the 
common envelope is greatly enhanced. The remaining 4 
systems (4U 1543-47, GRO J1655-40, GS 2023+338, and 
GRS 1915-1-105) have either intermediate-mass donor stars 
(i.e. 2 < Md < 5 Mq) or orbital periods longer than 2.5 days, 
thereby allowing for primordial secondaries of at least inter- 
mediate mass, and substantial mass loss or evolution of the 
secondary to its present status as the donor star. It is ou the 
evolution of these latter two categories, with particular em- 
phasis on GRS 1915-1-105, that we focus this work (for other 
recent discussions of intermediate-mass black-hole binaries 
see Kalogera 1999; Brown et al. 2000; LBW). 

In addition to the common-envelope ejection mecha- 
nism, another major uncertainty in the modelling of black- 
hole binaries is the initial mass of the black hole which is 
caused by uncertainties in the theory of both single and 
binary stellar evolution. Some of the key factors that de- 
termine the maximum initial black-hole mass are (1) the 
minimum initial mass above which a star leaves a black-hole 
remnant (mostly believed to be in the range of 20 - 25 Mq ; 
Maeder 1992; Woosley & Weaver 1995; Portegies Zwart et 
al. 1997; Ergma & Fedorova 1998; Ergma & van den Heuvel 
1998; Brown et al. 2000; Fryer & Kalogera 2001; Nelemans 
& van den Heuvel 2001; cf Romani 1992), (2) the minimum 
mass above which a single star loses its envelope in a stel- 
lar wind and becomes a helium/ Wolf- Rayet star, (3) the 
maximum radius of a single star before and after helium 
core burning, (4) the amount of wind mass loss in the Wolf- 
Rayet phase and (5) the fraction of the mass that is ejected 
when the black hole forms (for detailed recent discussions 
see Brown et al. 2000; Fryer & Kalogera 2001; Nelemans 
& van den Heuvel 2001). Generally one expects the most 
massive black holes to form from stars that have an initial 
mass close to the minimum mass above which a star loses 
its hydrogen-rich envelope in a stellar wind and becomes a 
Wolf-Rayet star, and where the common-envelope phase oc- 
curs near the end of the evolution of the massive primary 
(i.e. experiences case C mass transfer; Brown, Lee & Bethe 
1999; WcUstcin & Langcr 1999). This avoids a long phase 
where the mass of the helium star, the black-hole progeni- 
tor, is reduced by a powerful stellar wind, as typically seen 
from Wolf-Rayet stars, which would reduce the final helium- 



star mass and hence the maximum black-hole mass (see e.g. 
Woosley, Langcr & Weaver 1995)^. Unfortunately, the evo- 
lutionary tracks for massive post-main-sequence stars and in 
particular the maximum radius a star attains after helium 
core burning are rather uncertain (and generally inconsis- 
tent with observed distributions of stars in the Hertzsprung- 
Russell diagram; see e.g. Langer & Maeder 1995). Fryer & 
Kalogera (2001) have shown that initial black-hole masses as 
high as 15 M© can be obtained if either the parameter range 
for case C mass transfer is increased or the wind mass-loss 
rate in the helium-star phcise of the black-hole progenitor is 
reduced (see also Brown et al. 2001; Nelemans & van den 
Heuvel 2001; Belczynski & Bulik 2002; Pols & Dewi 2002) 

Only a few of the previous studios (e.g. LBW) have 
considered the possibility that the black hole may increase 
its mass substantially since its formation by mass transfer 
from the companion star. It is one of the purposes of this 
paper to demonstrate that accretion from a companion star 
can substantially increase the mass of a black hole and spin 
it up in the process and that the present mass may not 
be representative of the initial black-hole mass. A closely 
coupled result is that the observed donor star masses may 
be substantially lower than their initial mass. 

The paper is structured as follows. In § 2 we present de- 
tailed binary population synthesis calculations to show how 
the formation rate of black-hole binaries and the distribution 
of the secondary masses depend on the structure of massive 
supergiant envelopes and the modelling of common-envelope 
ejection. In § 3 we discuss the results of extensive binary evo- 
lution calculations for black-hole binaries with intermediate- 
/high-mass secondaries, which we then apply in § 4 to ob- 
served systems, in particular GRS 1915-1-105, ultraluminous 
X-ray sources and Cyg X-1. Finally, in § 4 we reexamine the 
standard formation scenario for low-mass black-hole binaries 
to understand why such systems appear to be so plentiful in 
Nature. 



2 BINARY POPULATION SYNTHESIS 
MODEL 



2.1 Assumptions of the Model 

In this work we consider only black-hole binaries that de- 
scend from a primordial binary pair of stars. Black-hole bi- 
naries that may form dynamically in globular clusters are 
left for another study. We start our investigation of black- 
hole binaries, which includes (initially) intermediate- and 
high-mass donors, by carrying out a Monte Carlo popula- 
tion study aimed at producing systems at the end of the 
common-envelope phase. This will provide guidance for the 
second part of our study where we follow in detail the X-ray 
binary phase which involves mass transfer from the donor 
star to the black hole. 

^ It should also be noted that in the formation of some blaek 
holes (e.g. the black hole in Nova Scorpii) a significant fraction of 
the mass of the helium star is ejected in the supernova explosion 
in which the black hole formed (Podsiadlowski et al. 2002; LBW). 
Thus the final helium-staj: mass strictly provides only an upper 
limit on the black-hole mass. 
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The assumptions and ingredients that we adopt for the 
population study are listed below. For prior studies of the 
formation and evolution of black hole binaries sec e.g. Ro- 
mani (1992), Portegies Zwart, Verbunt & Ergma (1997), 
Ergma & van den Heuvol (1998), Ergma & Fcdorova (1998), 
Kalogera (1999), Brown et al. (2000) and Fryer & Kalogera 
(2001); also see Kalogera & Webbink (1998) for a related, 
detailed study of the theoretical constraints on the formar- 
tion of neutron-star X-ray binaries. 

We utilize a simple power-law distribution for the initial 
mass function (IMF) for the primary stars in primordial bi- 
naries. Specifically, wc take dN/dMp oc Mp^ with x = 2.35 
(Salpctor 1955). Since wc consider a relatively small range 
of primary masses, the results are not very sensitive to the 
particular choice of x or to the fact that the IMF flattens 
significantly toward lower masses (sec e.g. Miller & Scalo 
1979). Only primordial binaries with mass Mp > 25 Mq are 
considered as progenitors of black holes. This lower limit is 
somewhat uncertain, but is consistent with current models 
of massive stars and the modelling of supernova explosions 
(e.g. Wooslcy & Weaver 1995; Fryer & Kalogera 2001). We 
also somewhat arbitrarily take the upper mass limit for the 
primary to be Mp < 45 Mq. 

The mass of the secondary star, Ms, in the primor- 
dial binary is chosen from a fiat mass ratio distribution, i.e. 
f{q) = 1, where q = Ms/Mp. There is a large degree of un- 
certainty in the actual distribution of mass ratios in binaries; 
however, our choice reflects the simple fact that a significant 
fraction of high-mass stars are observed to have high-mass 
companions (see e.g. Garmany, Conti & Massey 1980). At 
the low-mass end, only secondaries with mass Ms > 0.5 Mq 
are retained in the population synthesis. 

The orbital-period distribution of primordial binaries is 
taken to be constant in logPorb, where Porb is the orbital 
period (see e.g. Abt & Levy 1978). While orbital eccentric- 
ity among primordial binaries might be expected to have 
a distribution that increases linearly with e (Duqucnnoy & 
Mayor 1991), we have simply taken all orbits to be circular. 
One motivation for this choice is that the tidally circularized 
orbital radius of an eccentric binary is (1 — e^)a, which we 
note is simply linearly proportional to a, the initial orbital 
semimajor axis. Thus, the distribution in circularized orbital 
radii would be the same as that for the semimajor a^xes, re- 
gardless of the distribution in e, as long as that distribution 
is independent of a. 

The evolution of the primary star as it expands toward 
filling its Roche lobe is followed with the Hurley, Pols & Tout 

(2000, hereafter HPT) code. The evolution includes wind 
mass loss from the primary according to the prescription 
of Nieuwenhuijzen & de Jager (1990). These stars have fully 
developed core masses given by Mcore — 0.12 Mp '^^ (in solar 
masses; see e.g. HPT). Typical wind mass losses prior to 
the primary traversing the Hcrtzsprung gap (hereafter HG) 
amount to only ~ 2 — 6 Mq ; however, much larger losses can 
be sustained once stellar radii of > 1000 R© are reached. In 
this latter regard, we compute the orbital widening due to 
the stellar wind mass loss according to da/a = —dM^/M\y, 
where the fractional change in orbital separation is equal 
to the wind mass loss from the primary (dMw) divided by 
the total mass of the binary, Mb. This expression is based 
only on the assumption that the speciflc angular momentum 



carried away by the wind is equal to that of the primary star 
in the binary orbit (see, however, the discussion in § 4.4). 

An interesting possibility occurs if the initial Roche-lobe 
radius, 7?l, of the primary falls in the range Phg < -Rl < 
Pmax, where Rug is the stellar radius at the end of the 
Hertzsprung gap, and i?max is the ma^ximum radius that the 
primary can attain. It is in this phase of the evolution that 
the primary can lose a substantial fraction of its envelope 
mass in a stellar wind, thereby requiring less orbital bind- 
ing energy to eject the remaining envelope of the primary 
once mass transfer commences (see the discussion below). 
However, as the primary expands and loses mass in a wind, 
the mass loss also causes the orbit to expand (as described 
above). It it therefore not obvious whether the expanding 
star can catch up with its Roche lobe. To quantify this is- 
sue, wc first define a Roche-lobe index due to wind loss: 
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as well as a stellar index associated with wind mass loss: 

d\nR 



dlnMb/,, 



(0.023 + 0.00086 [30 - Mhg]) Mb(2) 



where Mb is the instantaneous total mass of the binary (in- 
cluding wind mass loss), Mhg is the mass of the primary 
at the end of the Hertzsprung gap (both in solar masses), 
and tl is the Roche-lobe radius of the primary in units of 
the orbital separation. In both cases, a minus sign has been 
subsumed into the definition of dMy, since it is always neg- 
ative. The above expressions for ^*,w and ^l,w were derived 
from fits to the stellar models in HPT and the expression 
for the Roche- lobe radius of Plavec (1968), respectively. If 
^*,w < Cl,w when the primary expands past the HG, then it 
will never catch up with its Roche lobe since ^,,w decreases 
faster with mass loss than does ^l.w (this catch-up problem 
has first been identified by Kalogera & Webbink [1998] in 
their study of the formation of neutron-star X-ray binaries). 
Wc have found that even if the reverse inequality holds when 
the primary expands past the HG, it is still extremely rare 
for the star to overtake its Roche lobe. 

If the primary does evolve to the point of overfiowing 
its Roche lobe (almost always before the end of the HG), we 
use a simple prescription for deciding whether a common- 
envelope phase ensues. If the primary has evolved at least to 
the start of the Hertzsprung gap and the mass ratio Mp/Ms 
exceeds 2.0, or if the primary is beyond the HG and the 
mass ratio exceeds a value of 1.2, we assume that a common 
envelope will occur. For other conditions we take the mass 
transfer to be quasi-conservative and stable. These latter 
systems are not very common and do not lead to the type 
of black-hole binary that we are considering. 

A commonly used prescription was employed to de- 
termine the orbital separation after the common-envelope 
phase (see e.g. de Kool 1990; Dewi & Tauris 2000): 
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where the subscripts / and i denote the final and initial 
values, respectively, Mc and Mo are the core and envelope 
masses of the primary, respectively, acE is the efficiency with 
which the orbital binding energy can liberate the envelope of 
the primary, is the energy of the envelope of the primary 
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Figure 1. The envelope structure parameter A as a function of stellar radius for different masses as indicated after hydrogen has been 
exhausted in the core. In the panels on the left, A only includes the gravitational binding energy, while on the right A includes both the 
gravitational and the thermal energy (similar to Dewi &; Tauris 2000). The dotted curves are calculated without inclusion of a stellar 
wind. Note that in this case A always becomes large at the largest radii attained by the models. 



in units of —GMpMc/Rp, and rL is the dimensionless Roche- 
lobe radius of the primary. 

Appropriate values of the A parameter are derived from 
stellar structure and evolution calculations (with similar 
assumptions as in HPT and wind mass loss according to 
Nieuwenhuijzen & de Jager 1990) that we have carried out 
specifically for this work. These results are shown in Fig- 
ure where the definition for A in the left panel includes 
only the gravitational binding energy (Agr), while in the 
panel on the right it also includes the thermal and the ioniza- 
tion energy (Atot; see Han, Podsiadlowski & Eggleton 1994; 
Dewi & Tauris 2000). We have compared these curves with 
those computed by Dewi & Tauris (2000, 2001) and found 
them to be in reasonable agreement at similar evolutionary 
stages, although the evolutionary tracks are slightly differ- 
ent. We note that the value of A itself depends to some degree 
on the definition of the core - envelope boundary (see Han 
et al. 1994 and Tauris & Dewi 2001). Here we have defined 
the core mass as the central mass that contains 1 Mq of hy- 
drogen (in contrast Dewi & Tauris [2000] defined it as the 
central mass which includes 10 per cent of the total mass 
of hydrogen). In agreement with Tauris & Dewi (2001) we 
find that in some evolutionary phases the value of A can 
be a rather sensitive function of the chosen core - envelope 
boundary; this is particularly true for very evolved super- 
giants near the end of their evolution. Unfortunately, this 
uncertainty in the definition of A cannot easily be resolved 
without a better understanding of the CE ejection process. 

Values of Atot that are appropriate to stars in the mass 
range 25 — 45M0 and the later phases of their evolution (i.e. 
in or beyond the HG) lie in the range of 0.01 ^ A < 0.06. Note 
that for a 20 M© model A increases to a value ~ 1 near the 
very end of the evolution. This occurs when the star ascends 
the asymptotic-giant branch and develops a deep convective 
envelope, creating a steep chemical gradient (/^-gradient) be- 
low the convective envelope (i.e. establishes a typical giant 
structure). The more massive stars experience a supernova 
before this phase is reached. Since these results depend on 
the assumptions in the stellar modelling, we also performed 
a series of calculations without any wind mass loss, shown 



as dotted curves in Figure The models now achieve sig- 
nificantly larger radii, and A increases dramatically near the 
end of the evolution in all cases. This demonstrates that the 
behaviour of A is also quite sensitive to uncertainties in the 
stellar modelling itself, adding another uncertainty to the 
problem. In view of these uncertainties we adopted a value 
of A in this work that is held constant throughout the evo- 
lution of a star. We were careful, however, to test the full 
range of plausible values for A (0.01 — 0.5). 

The value for acE was simply taken to be unity. Since 
the two parameters acE and A appear as a product in the ex- 
pression for the post-CE orbital separation, the uncertainty 
in one can, to some degree, be incorporated into the uncer- 
tainty in the other. A value of acE ~ 1 was motivated by a 
number of recent empirical studies of the efficiency of the CE 
ejection process which have suggested that the CE ejection 
process must be very efficient; see e.g. the recent study of 
sdB stars in compact binaries by Han et al. (2002a,b) (these 
are short-period binaries which have formed through a very 
well defined CE channel and are hence particularly suitable 
for studying the CE ejection process empirically). In their 
best model, these authors found acE ~ 0.75 and that ~ 75 
per cent of the thermal energy of the envelope had to be 
used in the ejection process to explain the observed orbital- 
period distribution of compact sdB binaries. However, this 
empirical study, as well as all other previous ones, strictly 
apply only to systems where the donor is a giant with a 
convective envelope but not to systems which experience a 
CE phase when the donor star is in the Hertzsprung gap 
and has a radiative envelope. Such stars have fundamen- 
tally different internal structures and are much more cen- 
trally concentrated than more evolved (super-)giants with 
convective envelopes (see e.g. Fig. 2 in Podsiadlowski 2001) 
(this is e.g. reflected in the low value of A for stars in the 
Hertzsprung gap in Fig. 1). In such systems, the core - enve- 
lope separation may not be distinct enough to allow envelope 
ejection even if enough energy is available in principle (see 
e.g. Taam & Sandquist [2000] and Kalogera [2002; private 
communication]); hence a value of acE ~ 1 may not be an 
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appropriate one for systems that experience a CE phase in 
the Hertzsprung gap. 

Immediately after the common envelope phase has oc- 
curred, we check whether the secondary star is overfilling 
its Roche lobe. If so, we assume that the secondary merges 
with the core and do not follow the binary further. 

The amount of wind mass loss from the exposed H- 
exhausted core of the primary as it evolves toward core 
collapse is quite uncertain. To cover a full range of possi- 
bilities we adopt the following somewhat ad hoc prescrip- 
tion for the mass of the black hole that eventually forms 
from the collapse: Mbh = 6M0 -|- (Mho - 6 M0)7^, where 
Mhb is the mass of the newly exposed H-exhausted core, 
and TZ is sl uniform random variable between and 1. 
With this prescription we produce black holes of mass over 
the range of ~ 6 — 20 Mq. The orbital separation at the 
end of this wind loss phase is taken to be a factor of 
(Mhg + Ms)/(AfBH + Ms) larger than the separation at the 
end of the common-envelope phase (see the expression above 
for the change in orbital separation with wind mass loss). We 
specifically chose this procedure to obtain black-hole masses 
that are reasonably consistent with the observed ones (e.g. 
Fryer & Kalogera 2001; van den Heuvel 2001; LBW), but 
also to avoid some of the theoretical uncertainties that lead 
to these masses (e.g. the wind mass-loss rate in the helium- 
star phase, the question of case B vs. case C mass transfer; 
see the discussion in § 1). If it is indeed necessary that sys- 
tems with relatively massive black holes have experienced 
case C mass transfer in the past (e.g. Brown et al. 1999; Well- 
stein & Langer 1999; Langer 2002 [private communication]), 
our procedure implicitly assumes that those evolutionary 
tracks for single stars which presently do not allow case C 
mass transfer for stars more massive than ~ 20 M0 are not 
correct (also see Fryer & Kalogera 2001). Note also that our 
procedure does not produce relatively low-mass black holes 
with masses < 6 Mq which could have a somewhat different 
evolution from the black-hole systems studied in this paper 
(see Fryer & Kalogera 2001 and Beer & Podsiadlowski 2002). 

When the evolving He star undergoes core collapse, we 
assume that the entire mass of the He star collapses into the 
newly formed black hole (but see LBW), and that there is 
no natal "kick" imparted to the black hole. Much has been 
inferred about natal kicks to neutron stars from the proper 
motions of radio pulsars and from studies of numerous in- 
dividual binary systems containing neutron stars (see e.g., 
Lyne & Lorimer 1994; Brandt & Podsiadlowski 1995; Ver- 
bunt & van den Heuvel 1995; Hansen & Phinney 1997; van 
den Heuvel et al. 2000; Arzoumanian, Chernoff & Cordes 
2002). Yet, relatively little is understood about the physical 
origin of these kicks or the conditions under which kicks are 
developed (see e.g. Janka & Miiller 1994; Fryer, Burrows & 
Benz 1997; Spruit & Phinney 1998; Fryer & Heger 2000; Lai 
2000; Pfahl et al. 2002). Unfortunately, even less is known 
about the nature of natal kicks delivered to forming black 
holes. The spatial distribution and the kinematics of the 
majority of black-hole binaries appear to be consistent with 
the assumption that no asymmetric kicks are imparted to 
the black hole (Brandt, Podsiadlowski & Sigurdsson 1995; 
White & van Paradijs 1996), although there is at least one 
clear exception (the black hole in Nova Sco; Brandt et al. 
1995; Fryer & Kalogera 2001; Podsiadlowski et al. 2002a; 



for a different view see Nelemans, Tauris & van den Heuvel 
1999). 

In order to compute an absolute formation rate for 
black-hole binaries in the binary population synthesis (here- 
after BPS; see next section), we take the birth rate of stars 
with mass > 8 Mq to be 1 per 100 years, to match the Galac- 
tic supernova rate (of non-type la SNe; see e.g. Cappellaro, 
Evans & Turatto 1999). This rate, coupled with the assumed 
slope of the IMF, yields a maximum potential formation rate 
for black hole primordial binaries of ~ 6 x 10~* yr~^; this 
includes an approximate fraction for stars born in binaries 
of ~ 1/2. 

Finally, we note that once the incipient black-hole bi- 
naries have been produced, the luminous X-ray phase will 
not start until the secondary has evolved to either produce 
a strong stellar wind or overflow its Roche lobe. In the BPS 
portion of this work we carry the calculations only to the 
point of the successful ejection of the common envelope or 
merger of the secondary with the core of the primary. We 
also apply a very simplistic check for the dynamical sta- 
bility of mass transfer when the donor star ultimately com- 
mences mass transfer via Roche-lobe overflow onto the black 
hole. Here we required that Ma ^ Mbh; however, we note 
that in the detailed binary evolution calculations presented 
in §3, the stability of the mass transfer is calculated self- 
consistently in each evolution step, and that even higher 
mass donor stars may result in stable mass transfer. The 
successful incipient black-hole binaries are tabulated and 
statistical results are extracted. 

Many of the above assumptions regarding the forma- 
tion of black-hole binaries are discussed in further detail by 
Romani (1992), Portegies Zwart, Verbunt & Ergma (1997), 
Ergma & van den Heuvel (1998), Ergma & Fedorova (1998), 
Kalogera (1999), Brown et al. (2000), and Fryer & Kalogera 
(2001). One substantive difference in the assumptions that 
we make as compared to those utilized in prior work is that 
we do not set any a priori restrictive upper limit on the mass 
of the donor star in the incipient black-hole binary. Also, the 
values for A that we use are considerably smaller (by up to 
a factor of ~ 20) than the values employed in most earlier 
work. Smaller values of A require the expenditure of more 
orbital binding energy to successfully eject the common en- 
velope. Finally, we do not require that the donor star be 
unevolved when it commences mass transfer onto the black 
hole (cf. Kalogera 1999). This is due to the fact that we 
explicitly include wide binaries containing black holes. 

2.2 Binary Population Synthesis Results 

Each binary population synthesis run (BPS) was typically 
started with ~ 10*^ primordial binaries, all of which have 
Mp > 25 M0 . Each binary is followed using the prescrip- 
tions, assumptions, and algorithms specified in the previous 
section, until it either becomes an incipient black-hole X-ray 
binary or goes down an alternate evolutionary path. The 
relevant binary parameters are stored for each "successful" 
system. 

Illustrative results from a BPS run with the A- 
parameter set equal to 0.5, a conventionally used value, are 
shown in Figure Each dot in the Porb — Md,i plane rep- 
resents a single incipient black-hole X-ray binary just after 
the black hole has been formed; Porb is the orbital period 
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Figure 2. Distribution of initial orbital period and initial secondary mass for black-hole binaries at the beginning of mass transfer for 
different assumptions about the common-envelope (CE) structure and ejection efficiency, characterized by the parameters A and acE, 
respectively, where = GMeMp/XRp is the binding energy of the envelope. The value of ocE wa^s fixed as 1, while two different 

illustrative values of A were used to produce the right and left panels. The solid circles indicate the initial positions of our detailed binary 
evolution sequences which are discussed in §3. 
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Figure 3. Formation rate of black-hole binaries as a function of the envelope structure parameter A. The different curves assume a 
different maximum mass for the secondary (in M©), as indicated along the respective curves. 



and Md,i is the initial mass of the donor star. The number 
of primordial binaries used in this run was limited to only 
1.5 X 10® in order to keep the density of points legible in 
the figure. The sharp lower boundary in Figure results 
from the fact that systems with shorter orbital periods have 
merged, i.e. the donor star overflowed its critical potential 
lobe at the end of the common-envelope phase. The filled 
circles are starting models for the detailed evolutionary cal- 
culations of the binary X-ray phase presented in § 3. 



In all, there are about 25,000 dots in Figure each 
representing the successful formation of an incipient black- 
hole binary X-ray source. The expression we use to convert 
this number, A'^, to a formation rate is as follows: 

where A^o is the starting number of primordial binaries (with 
Mp > 25 Mq), the factor of 1/2 is an approximation to the 
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fraction of stars born in binaries, the factor of 1/8 represents 
the fraction of stars capable of collapsing to a neutron star 
or black hole that comes from the mass range 25 — 45 M©, 
and is the Galactic supernova rate (of non-type la SNe; 
see discussion above). The factor of 1/8 depends somewhat, 
but not sensitively, on the slope of the IMF. The particular 
example shown in Figure then implies a formation rate 
for such systems of ~ 1 x 10~^ yr~^- 

A second illustrative example of a BPS run is shown 
in Figure |5j3 for the case where the A-parameter has been 
reduced to 0.08. The number of primordial binaries in this 
run was 10 times that used in Figure |2^. In spite of this, 
the number of systems represented in Figure is actually 
fewer than in Figure |2t due to the highly reduced formation 
efficiency associated with the lower value of A. Note that, in 
addition to the overall drop in formation efficiency, there is 
a complete dearth of systems with donor masses of < 5 Mq . 
This results from the fact that for small values of A (i.e. 
more envelope binding energy) the systems with lower-mass 
donor stars do not have sufficient orbital energy to unbind 
the common envelope. 

We have systematically investigated the effect of A on 
the formation efficiency of incipient black-hole binaries. The 
BPS code was run for 50 values of A in the range of 0.01 — 0.5, 
in steps of 0.01, and a formation efficiency computed for 
each. The results, shown in Figure |5] are further broken 
down according to the donor masses contributing to the pop- 
ulation. The top curve shows the formation rate for systems 
with all donor masses (< 15 M©) computed according to 
equation (4). Note the abrupt drop in formation rate for 
values of A<0.1. The sequence of curves below this is for 
the formation rate of black-hole X-ray binaries with donor 
masses < 8, < 5, < 3, < 2, and < IMq, respectively. The 
small shoulder on the curve for < 15 Mq at very small val- 
ues of A results from Roche-lobe overflow when the primary 
has evolved beyond the HG and has lost a significant amount 
of its envelope to a stellar wind. Two immediate conclusions 
to be drawn from the results in Figure are: (1) black-hole 
binaries with low incipient donor masses (i.e. < 1 Mq) have 
extremely low formation rates unless either seemingly un- 
physical (large) values of A are invoked or much lower-mass 
primary stars are able to form black holes (see e.g. Porte- 
gies Zwart et al. 1997; Kalogera 1999); and (2) the forma- 
tion rates for black- hole binaries with incipient intermediate- 
mass donor stars can be quite substantial provided that the 
A-parameter is not smaller than ~ 0.1. 

In addition to the BPS results shown in Figure El we 
have carried out some additional studies of how the the- 
oretical uncertainties in the values of A affect our conclu- 
sions. As can be seen from Figure the values of A depend 
systematically on the stellar mass, with the larger values 
(smaller envelope binding energies) generally being associ- 
ated with the lower masses. Since, to this point, we have 
excluded primary masses below 25 Mq , and these have the 
largest values of A, we have carried out an additional se- 
quence of BPS calculations with the inclusion of 20 - 25 Mq 
primaries to test whether these would significantly enhance 
the formation rate of black-hole binaries. Here, we simply 
recomputed the formation rates, as in Figure |3| but now 
with the minimum primary mass set to 20 Mq. Again, for 
the purposes of the BPS calculations, we took the value of 
A to be fixed for primaries of all masses and evolutionary 



states, and studied how the formation rates changed with a 
systematic variation in the adopted value of A. The result is 
that for any given value of A, the rates are larger than those 
shown in Figure |3 by a factor typically limited to ~ 1.5. 
This is a bit less than might be expected simply from the 
increased primary mass range - weighted by the Salpeter 
mass function. The smaller than anticipated increase is due 
to the fact that many of the systems produced by ~ 20 Mq 
primaries have relatively low-mass black holes and the sub- 
sequent mass transfer is unstable. 

Finally, we constructed from the results of Figure a 
simple dependence of the A parameter on primary mass: 
A = Ao exp(— 0.074 [Afp — 20 Mq]), where Ao is a constant 
to be supplied. For the 'plateau' regions in Figure the 
value of Ao would be approximately 0.08. We then produced 
a sequence of BPS runs, including primary masses down to 
20 Mq , where the value of Ao was varied systematically over 
the range of 0.01-0.5. Thus, for any given choice of Aq, in 
a given BPS run, the values of A used scale with primary 
mass as in the exponential expression given above. The net 
result of this study, as expected, is that the formation rates 
are all systematically reduced compared to the values shown 
in Figure |3 where a fixed value of A was utilized in each 
BPS run. The reason for this reduction in formation rates 
is that now Aq represents an upper limit to the value of A 
(corresponding to 20 Mq stars), with all higher-mass stars 
having smaller values of A. 

The conclusion drawn from these experiments, that ex- 
tend the primary masses down to 20 Mq and compute A 
from a prescription that depends on the primary mass, is 
that the rates shown in Figure |3 are, to within a factor of a 
few, likely to be upper limits. 



3 BINARY EVOLUTION CALCULATIONS 

In this section we present a series of binary stellar evolution 
calculations to illustrate the evolution of black-hole binaries 
with intermediate- and high-mass secondaries. In § 3.1 we 
briefiy describe the binary evolution code and the main as- 
sumptions used in the calculations, and in § 3.2 we present 
the main results of the calculations. In the subsequent sec- 
tion we will apply these results to a variety of systems. 

3.1 Description of the code and principal 
assumptions 

All our calculations were performed with a standard 
Henyey-type stellar evolution code (Kippenhahn, Weigert 
& Hofmeister 1967) which we have used in various similar 
investigations in the past and which is described in detail in 
Podsiadlowski, Rappaport & Pfahl (2002b; hereafter PRP). 
It uses up-to-date stellar OPAL opacities (Rogers & Iglesias 
1992), complemented with those by Alexander & Ferguson 
(1994) at low temperatures. In all calculations we assumed 
an initial solar composition with X = 0.70 and Z = 0.02, 
took a mixing-length parameter a — 2.0, and included 0.25 
pressure scale heights of convective overshooting from the 
core (Schroder, Pols & Eggleton 1997; Pols et al. 1997). 

Since the accreting object is a black hole, we need to 
follow the change in the accretion efficiency and the change 
in the spin angular momentum of the black hole as its mass 
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Figure 4. Selected properties of tlie evolutionary sequences for black-hole binaries with initially unevolved secondaries of 2 to 17 Mq 
(see Table 1) as a function of time since the beginning of mass transfer (with arbitrary offset). The initial mass of the black hole is 10 
in all sequences. Top left panel: black-hole spin parameter (J/A/^); top right panel: orbital period; bottom left panel: black-hole mass 
(dashed curves) and secondary mass (solid curves); bottom right panel: mass-transfer rate. The shaded regions in each panel indicate 
the period range of 30 to 40 d (similar to the orbital period of GRS 1915+105 with Porb = 33.5 d). 



increases. In our calculations we assume for simplicity that 
the black hole is initially non-rotating. If we further assume 
that the efficiency, rj, with which the black hole radiates is 
determined by the last stable particle orbit, then the black- 
hole luminosity can be written as 

i = 77MaccC^ (5) 

where Mace is the black-hole mass-accretion rate (as mea- 
sured by an observer at infinity), c is the speed of hght and 
?7 is approximately given by 




for Mbh < \/6A/bh- Here, M^jj and Mbh are the initial 
and the present gravitating mass-energy of the black hole, 
respectively (Bardeen 1970; see also King & Kolb 1999). 
In none of our evolutionary sequences does Mbh exceed 
\/6M£jj. Over this interval rj ranges from ~ 0.06 — 0.42 (but 
see footnote 2). If the black hole were born with significant 



The Formation and Evolution of Black-Hole Binaries 9 



rotation (as argued e.g. by LBW), all of these expressions 
would have to be modified accordingly. 

This luminosity needs to be compared to the Eddington 
luminosity at which the radiation pressure force balances 
gravity 

AtvGMbhc 

-t^cdd ~ , (I) 

K 

where G is the gravitational constant and k is the opac- 
ity assumed to be due to pure electron scattering, i.e. 
K — 0.2(1 + X)cm^g~^ for a composition with hydro- 
gen mass fraction X (e.g. Kippenhahn & Weigert 1990). 
Equating L^dd to L in equation (|SJ then defines the Edding- 
ton mass-accretion rate, i.e. the maximum accretion rate at 
which gravity can overcome radiation pressure (for spherical 
accretion): 



Medd 



47rGA/BH 
Kcrj 



(8) 



2.6 X lO~''M0yr'^ 



Mb 



10 Mf 



(9) 



In most calculations we assume that any mass transferred in 
excess of the Eddington accretion rate is lost from the sys- 
tem, carrying with it the same specific angular momentum 
as the orbiting black hole, while the rest of the mass, reduced 
by the fractional rest mass energy lost in the radiation, is 
accreted by the black hole. 

As the black hole accretes mass and angular momen- 
tum, its spin parameter, a = J/M^, increases according to 



(I) 



1/2 tfO 



Mbh 



1/2 



(10) 



for Mbh < V6A/£h (see e.g. Bardeen 1970; Thorne 1974; 
King & Kolb 1999) ^ 

Our calculations include, as a default, orbital angular 
momentum losses via magnetic braking if the secondary has 
a convective envelope (see e.g. Verbunt & Zwaan 1981; Rap- 
paport, Verbunt & Joss 1983), and by gravitational radia- 
tion; these affects are only important in some of the calcu- 
lations with relatively low-mass donor stars in § 3.2 (also 
see § 4.4). In most cases mass transfer occurs either on a 
thermal timescale or is driven by the nuclear evolution of 
the secondary. 



and 13 and in the top part of Table 1. All calculations were 
terminated either when the secondary became degenerate 
or when it became detached during helium core burning. In 
some cases, the secondary would have filled its Roche lobe 
again after helium core burning and ascended the asymp- 
totic giant branch. In this relatively short-lived phase the 
properties of the systems would be similar to those on the 
ascent of the first giant branch and the orbital period would 
continue to increase somewhat (but typically by less than a 
factor of 2). In all cases with initially unevolved secondaries, 
the secondaries end their evolution as white dwarfs rather 
than in a supernova. 

The main behaviour of these binary sequences is not 
difficult to understand; it is mainly determined by the ini- 
tial mass ratio of the components (g = Aid /Mbh', see e.g. 
Ritter 1996; Kalogera & Webbink 1996; PRP for detailed re- 
cent discussions) . Generally the mass-transfer rates (bottom 
right panel in Fig. 01 at any point in the evolution are higher 
for higher initial secondary masses (mass ratios). For rela- 
tively low-mass secondaries, the secondaries always remain 
close to thermal equilibrium and mass transfer is entirely 
driven by the nuclear evolution of the secondary. The mass- 
transfer rate tends to dip near the end of the secondary's 
main-sequence phase, where the secondary may even become 
detached temporarily, and then increases sharply as the sec- 
ondary ascends the giant branch where the evolution is de- 
termined by the rate at which the hydrogen-burning shell 
advances through the star. For the more massive secondaries 
{q > 1), the mass-transfer timescale becomes shorter than 
the thermal timescale of the secondary and mass transfer 
occurs initially on the thermal timescale of the secondary's 
envelope. For the most massive secondaries this leads to a 
sharp initial spike in the mass-transfer rate, reaching a peak 
of M ~ 10~* M0 yr~^. Note that after this initial high mass- 
transfer rate, the evolution of all sequences with M > 12 M0 
becomes essentially uniform and independent of the initial 
mass of the secondary^. 

For secondaries with initial masses larger than ~ 20 Mq 
(not shown in the figures; but see § 4.3), the initial mass- 
transfer rate would become so high that mass transfer would 
become unstable (i.e. experience a delayed dynamical insta- 
bility; Hjellming & Webbink 1987; PRP) and the secondary 
would then be likely to engulf the black hole leading to a sec- 
ond common-envelope phase and the spiral-in of the black 
hole in this envelope''. The further evolution in this case is 



3.2 Results of the Binary Evolution Calculations 

To illustrate the evolution of black-hole binaries we first per- 
formed a sequence of models where we varied the mass of 
the secondary from 2 to 17 Mq. In these calculations the 
secondary was an initially unevolved (i.e. zero-age main- 
sequence) star and the black-hole had an initial mass of 
10 Mq. Furthermore, we assumed that the mass-accretion 
rate onto the black hole was Eddington-limited (eq.|5J. Some 
of the main results of these sequences are shown in Figures^] 

^ For a maximally spinning Kerr black hole a nominally ap- 
proaches unity, but is probably limited by a counteracting torque 
due to disk radiation swallowed by the black hole to be ~ 0.998 
(Thorne 1974). It also limits the efficiency rj in equation JHJ to 
~ 0.30. 



^ This is only true for secondaries that are initially relatively un- 
evolved which after the initial high mass-transfer phase behave 
like lower-mass secondaries. If the secondaries have already es- 
tablished a very non-uniform chemical composition profile, the 
subsequent evolution can be drastically altered (see e.g. Podsiad- 
lowski &; Rappaport 2000). 

■* In the calculation with an initially unevolved 20 Mq star, the 
secondary never overfilled its Roche lobe by more than 14 per 
cent (for a short period of time). It is not entirely clear whether 
this necessarily leads to a spiral-in phase, in particular for stars 
with radiative envelopes (Podsiadlowski 2001). If it does not, and 
the system survives the initial phase of high mass transfer, the 
secondary becomes detached after it has reestablished thermal 
equilibrium (similar to Fig. 7 in PRP); but the secondary now 
has a mass of only 4 Mq and the subsequent evolution will mimic 
the evolution of a system with a 4 Mq secondary. It starts to fill 
its Roche lobe shortly after the hydrogen core-burning phase (i.e. 
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Figure 5. Shaded contours of black-hole spin parameter (J/M^) (top panels) and potential X-ray luminosity (bottom panels) (i.e. the 
luminosity assuming that all the mass transferred is accreted by the black hole) in the Porb^Mj (orbital period - secondary mass) 
plane (left panels) and in the H-R diagram (right panels). Shading for the black-hole spin (from light to dark): 0-0.25; 0.25-0.50; 
0.50-0.75; 0.75 — 1.00. Shading for the potential X-ray luminosity (from light to dark): transient X-ray sources; steady sources with Lx < 
1039 ergs /s); 10^^ < Lx < lO"*" ergs /s; Lx > lO"" ergs /s. The solid curves show the evolutionary tracks for the initially unevolved 
secondaries in Table 1 and Figure|3]with the initial masses as indicated. The dashed curves are contours of constant black-hole mass (10, 
11, 12, 13, 14, 15, 16 Mq from left to right [right panels], bottom to top [left panels]). The white strips indicate the position of systems 
with orbital periods between 30 and 40 d, i.e. have an orbital period similar to GRS 1915-1-105. 



presently rather uncertain; the black hole will most likely 
settle at the center of the secondary destroying its core or 
merging with it. The ultimate product will be a black hole 
with a possibly much larger mass than the initial mass if it is 



experience early case B mass transfer) and has parameters at an 
orbital period of 33.5 d consistent with those of GRS 1915-1-105 
(see § 4.1). 



able to accrete a substantial fraction of the secondary (e.g. if 
accretion occurs in the super-critical regime where radiation 
can be trapped in the flow; Houk & Chevalier 1991; Cheva- 
lier 1993). In addition, the black-hole may be surrounded by 
planet-mass objects or one or more low-mass stars that are 
likely to form by gravitational instabilities in the centrifu- 
gally supported disk left over from the collapsing envelope of 
the secondary (for more details of such a scenario see Pod- 
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Table 1. Selected Properties of Binary Evolution Sequences 
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2 


10 


0.56 


_ 


_ 


_ 


_ 


_ 


_ 


_ 


0.2.50 


11.63 


14 


3.1E+09 


3.1E+09 


_ 


_ 
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3 


10 


0.62 


0.573 


12.25 


4.2E-08 
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0.36 


0.337 
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8.6E+08 


8.3E+08 
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4 


10 


0.68 


0.802 


12.95 


1.2E-07 


3.644 


2.085 
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0.66 


0.439 


13.26 
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4.2E+08 


1.6E+08 


2.6E+08 


- 


- 


5 


10 


0.73 


0.995 


13.66 


2.1E-07 


3.666 


2.228 


0.441 


0.74 


0.556 


14.00 
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0.626 
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14.50 
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0.81 
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196 
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8 


10 


0.83 
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0.84 


0.783 


15.22 
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8.2E+07 


4.3E+05 
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10 
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2.8E-06 
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10 
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0.92 
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16.69 


134 


4.5E+07 


6.4E+05 


2.7E+07 


1.7E+07 


4.1E+05 


13 


10 


0.94 
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0.93 
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16.82 


124 
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1.21 
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4.5E-06 
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3.106 
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0.82 
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162 


4.0E+07 


1.5E+06 


3.7E+07 


1.3E+06 


6.1E+05 


8 


10 


1.86 


2.013 


13.77 


9.8E-06 


3.923 


3.463 


0.288 


0.75 


1.186 


13.84 


152 
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7.2E+05 


6.2E+06 


9.6E+06 


5.0E+05 


8 


10 


3.26 


2.634 


11.00 


3.7E-06 


3.936 


3.593 


0.528 


0.31 


1.564 


11.12 


138 


4.3E+06 


5.6E+05 


1.5E+06 


1.3E+06 


9.6E+05 


8 


10 


5.26 


3.268 


10.01 


2.6E-04 


3.850 


3.315 


0.700 


0.00 
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10.05 


131 
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9.7E+04 
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8 


10 


0.83 


1.162 


16.13 


1.3E-06 


3.762 


2.672 


0.331 


0.90 


0.729 


16.50 


106 


9.1E+07 


2.6E+06 


8.7E+07 


5.1E+05 


8.9E+05 


15 


10 


0.97 


1.383 


21.39 


3.3E-06 


3.941 


3.442 


0.238 


1.00 


1.163 


21.55 


48 


4.2E+07 


6.2E+05 


1.4E+06 


3.9E+07 


5.4E+05 


Note. - 


- Initial Parameters. M°, M 


jjj! initial 


masses 


of tliG secondary and tlie blaclt liole, P^^^^ 


initial orbital period. Parameters for GRS 1915 



(i.e. at Porb = 33.5 d). Mj, Mbii: masses of the secondary (the donor star) and the black hole, M: mass-transfer rate, Tes, L, A's„rt: effective 



temperature, luminosity, surface hydrogen mass fraction of the secondary, aspi„: black-hole spin parameter. Final parameters of the sequences. 
M\, Mgjj: final masses of the secondary and the black hole, P^b- fii^^l orbital period. X-ray lifetimes, tx-ray^ total lifetime of all X-ray phases, 
itrans: Overall lifetime of transient X-ray phases, t[<39], ([39,40], t\yia]'. lifetimes of steady X-ray phases with potential X-ray luminosities with 
Lx < 10^''ergs-\ lO^^ergs-' < ix < 10"° and Lx > iC^ergs^^ respectively 



siadlowski et al. 1995). If the black hole is indeed orbited by 
a low-mass stellar companion, the system may later appear 
again as an X-ray binary, but would now be classified as 
a low-mass black-hole binary. Note, however, that the sec- 
ondary should be chemically anomalous for its mass since it 
is likely to contain matter that underwent nuclear processing 
in a massive star (i.e. show evidence for CNO processing). 

Returning to Figure^] the left bottom panel shows that 
the black-hole mass (dashed curves) can increase substan- 
tially in these sequences, by up to ~ 7 M© (also see Table 1), 
even though we assumed that mass accretion onto the black 
hole was Eddington-limited. Most of this accretion takes 
place when the secondary is still burning hydrogen in the 
core since this phase lasts much longer than the subsequent 
giant phase and since the mass-transfer rates on the giant 
branch are generally much higher, often significantly super- 
Eddington. This has the consequence that a larger fraction 
of the transferred mass is lost from the system. 

Our finding that the black hole can accrete a fairly sig- 
nificant amount of mass is in apparent conflict with the re- 
sults of a related study by King & Kolb (1999) who found 
a much more moderate possible increase of the black-hole 
mass (<2.5M0) by estimating the maximum amount that 
can be accreted as the product of the Eddington accretion 
rate and the evolutionary timescale of the secondary. The 
differences in these results can be attributed to two factors. 
First King & Kolb (1999) used a characteristic value for the 
Eddington accretion rate (10~^ Mq yr~^ for a 10 Mq black 
hole) that is substantially lower (up to a factor of ~ 4) 
than the value given by equation JUJ in the early phase of 
the evolution when the black-hole radiation efficiency is low 
(7; ~ 0.06). A second factor is that, in particular for the 
more massive secondaries, the characteristic evolutionary 
timescale increases as the mass of the secondary decreases 
and as the secondary behaves like a less massive star. This 



increases the evolutionary timescale by between a factor of 
~ 2.5 for the least massive secondaries to ~ 4 for the more 
massive secondaries in our calculations (a factor of 100 in 
the 20 M© calculation). 

However, all of the calculations presented so far assume 
that the secondary is initially unevolved. As the results in § 2 
show, it is much more likely that the secondary was at least 
somewhat evolved at the beginning of mass transfer. Since 
this shortens the remaining evolutionary time in the hydro- 
gen core-burning phase, it reduces the amount of matter that 
can be accreted by the black hole. To illustrate this we have 
performed four evolutionary sequences for a secondary with 
an initial mass of 8Mq at different evolutionary stages (see 
'Evolved Sequences' in Table 1). In these calculations the 
secondaries had a hydrogen mass fraction in the core at the 
beginning of mass transfer oi X = 0.50, 0.30, 0.10 in the first 
three sequences, respectively, while in the fourth sequence 
the secondary already had developed a hydrogen-exhausted 
core of 0.7 M0 (i.e. experienced so-called early case B mass 
transfer). As expected, the final black-hole mass decreases 
systematically from 15.2 M© for the initially unevolved sec- 
ondary to only 11.1 Mq for the secondary near the end of the 
main-sequence phase. In the calculation where mass trans- 
fer starts after the main-sequence phase of the secondary, 
the black hole accretes only ~ 0.05 Mq. Nevertheless, even 
for an initially moderately evolved secondary, the black-hole 
mass can still increase quite substantially (by ~ 4M0). 

As the black hole accretes matter from the last stable 
orbit, it also accretes angular momentum and is spun up in 
the process. The top left panel in Figure |1] shows the time 
evolution of the black-hole spin parameter a (eg. IIOII . Even 
if the black hole was completely non-rotating initially (as we 
assumed in our sequences) and the accretion rate is Edding- 
ton limited, the black hole can be spun up substantially to 
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a spin parameter a ~ 0.9, where for a maximally rotating 
Kerr black hole a = 0.998 (Thorne 1974). 

In Figurel^we show the distribution (indicated by shad- 
ing) of the spin parameter (top panels) and the potential X- 
ray luminosity, defined below, (bottom panels) for selected 
evolutionary sequences of the secondary both in the orbital 
period - secondary mass plane (left panels), and in the H-R 
diagram (right panels). In these figures, the solid curves 
represent selected evolutionary sequences with initially un- 
evolved secondaries (as indicated). Note in particular how 
the evolutionary tracks for the most massive secondaries all 
converge and how the effective temperature in the red-giant 
phase increases systematically with initial secondary mass. 
The latter is in part caused by the much lower hydrogen 
abundance in the envelope of the initially more massive sec- 
ondaries (see Table 1). The dashed curves in the panels in- 
dicate the black- hole mass (from 10 to 16 M©). 

In the bottom panels of Figure|S|we have defined a ^po- 
tential X-ray luminosity' as the accretion luminosity from 
the black hole assuming that all the matter that is trans- 
ferred from the secondary is accreted by the black hole ra- 
diating at the appropriate efficiency r] (eq. (SJ, i.e. assuming 
that accretion is not Eddington limited (see e.g. Begelman 
2002). Because of the high mass-transfer rate, in particu- 
lar for the more massive systems, this potential X-ray lu- 
minosity can be as high ~ 10*^ ergs s~^ and systems can 
spend a large fraction of their X-ray-active lifetime at these 
high mass-transfer rates (see the X-ray lifetimes in Table 1 
and the further discussion in § 4.2). Figure |S| also shows 
where systems would be expected to be black-hole tran- 
sients (light shading). To decide whether a system exhibits 
transient behaviour, we utilized an expression very similar 
to equation (A5) of Vrtilek et al. (1990) for determining 
the outer disk temperature at r^j. If T^{rd) is found to be 
> 6500 K, we take the disk to be ionized, and therefore not 
subject to the standard thermal-ionization disk instability 
(Cannizzo & Wheeler 1984; van Paradijs 1996; King, Kolb 
& Szuskiewicz 1997; Lasota 2002). As expected from the 
behaviour of the mass-transfer rates (Fig. 2J, systems with 
relatively low-mass secondaries tend to be transient systems, 
and the potential X-ray luminosities increase systematically 
with initial secondary mass. 

When interpreting these results, several caveats are in 
order. The actual binary evolution calculations leading to 
the results discussed above assumed that any mass trans- 
ferred in excess of the Eddington rate was lost from the 
system; this, in turn, somewhat affects the evolution of 
the orbit and the mass-transfer rate itself. In order to es- 
timate how sensitive our results for the potential luminosity 
are to the assumption of Eddington-limited accretion, we 
performed two sequences where we assumed that accretion 
onto the black hole was not Eddington limited for secon- 
daries with an initial mass of 8 and 15 M© , respectively (see 
'Non- Eddington limited sequences' in Table 1). The 8M0 
sequence is only moderately affected since the mass-transfer 
rate is sub-Eddington for most of the evolution, while in 
the 15 M0 sequence the black-hole mass grows significantly 
larger, as expected, than in the Eddington-limited case. But 
even for the more massive secondary, the mass-transfer rates 
are typically within a factor of two at the same orbital pe- 
riod for the Eddington-limited and the non-limited case. 
We therefore conclude that the inferred potential X-ray lu- 



minosities are correct as computed to within a factor of a 
few. Since the Eddington luminosity in our systems with the 
highest mass black holes (~ I7M0) is 4 x 10^^ ergss"^, an 
observed X-ray luminosity as high as 10'"' ergss"^ would re- 
quire only a modest super-Eddington mass accretion rate of 
a factor of a few Afcdd . Even for systems with potential X-ray 
luminosities as high as 10*^ ergss"^, mass accretion has to 
exceed the Eddington accretion rate by typically less than a 
factor ~ 20. These super-Eddington accretion rates may be 
significantly reduced if beaming of the X-ray flux is impor- 
tant in these systems (King et al. 2001). Begelman (2002) 
has recently reexamined the problem of super-Eddington ac- 
cretion and concluded that in radiation-pressure dominated 
accretion discs super-Eddington accretion rates of a factor of 
10 to 100 can be achieved, due to the existence of a photon- 
bubble instability in magnetically constrained plasmas. In 
this context, it is worth pointing out that a number of X- 
ray binaries containing neutron stars are known to radiate 
substantially above the Eddington limit (SMC X-1, Levine 
et al. 1993; LMC X-4, Levine et al. 1991) by factors of up 
to ~ 5. It is believed that this is a consequence of the fact 
that the accretion onto the poles of the neutron star is fun- 
neled through a strong magnetic field, a process that is not 
directly applicable to black-hole systems. Nevertheless, it 
demonstrates that super-Eddington X-ray binaries exist in 
Nature (also see § 4.1). 

It is also worth pointing out that the mass-transfer rates 
obtained from our calculations are secular mass-transfer 
rates, i.e. represent an average over timescales much longer 
than the lifetime of X-ray astronomy. It is quite plau- 
sible that the mass-transfer rates fluctuate substantially 
about the secular mean even in systems that are not con- 
sidered 'transients' according to the disc- instability model. 
GRS 1915-1-105 may present an example for this. It became 
an X-ray source in 1992 (Castro-Tirado, Brandt & Lund 
1992) and has been a relatively steady source ever since 
(Sazonov et al. 1994; Greiner, Morgan & Remillard 1997; 

:/ /xte.mit.edu| |. Its behaviour is very different from 
the normal behaviour of soft X-ray transients and the sys- 
tem could probably be better classifled as a semi-persistent 
source. 



4 APPLICATIONS 
4.1 GRS 1915+105 

One of the best-studied black-hole binaries in the Galaxy 
is the microquasar GRS 19154-105 (see e.g. Castro-Tirado, 
Brandt & Lund 1992; Greiner, Morgan & RemiUard 1996; 
Mirabel & Rodriguez 1994). Greiner, Cuby & McCaughrean 
(2001) recently determined the orbital period of the system 
as 33.5 ± 1.5 d and obtained the black-hole mass function 
(/(M) = 9.5 ± 3.0 Mq). Based on their analysis, they find 
a black-hole mass of 14 ± 4 M©, substantially more massive 
than the masses inferred for the majority of black-hole tran- 
sients (see e.g. Table 1 in LBW and references therein; Orosz 
et al. 2002). In Figures|l|and|Sl we indicated the period range 
of 30 to 40 d, i.e. close to the period of GRS 1915-1-105, in 
all panels and in Table 1 we give some of the key system pa- 
rameters of all evolutionary sequences at an orbital period of 
33.5 d. The sequence with an initially unevolved secondary 
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Figure 6. X-ray luminosity, assuming Eddington-limited accretion, (left) and potential X-ray luminosity, assuming non-Eddington- 
limited accretion, (right) for the binary sequences with unevolved secondaries (see Fig. |4j as a function of time since the beginning of 
mass transfer. The shaded regions in each panel indicate the period range of 30 to 40 d (similar to the orbital period of GRS 1915-1-105 
with Porb = 33.5 d). 



of 2 M0 never reaches an orbital period of 33.5 d. The reason 
is that in this case (unlike the other cases), magnetic brak- 
ing becomes important, indeed the dominant mass-transfer 
driving mechanism, at the end of the donor's main-sequence 
phase; this causes a peak in M and a temporary shrinking of 
the orbit (see Fig. 4). We therefore added another sequence 
with an initial secondary of 2 Mq which has already evolved 
off the main sequence at the beginning of mass transfer (i.e. 
experienced case B mass transfer). Since the evolution of 
low-mass giants is entirely determined by the evolution of 
the core mass (and not the total mass), this sequence can 
also be considered representative for systems with low-mass 
secondaries in general (as assumed e.g. by Belczynski & Bu- 
lik 2002; Vilhu 2002 in their modelling of GRS 19154-105). 

As Table 1 shows, the masses of both components and 
the mass-transfer rate at an orbital period of 33.5 d in- 
crease systematically with the mass of the secondary, reach- 
ing a maximum for an initial secondary of ~ ISM©. The 
same applies to the effective temperature and the luminos- 
ity. Greiner et al. (2001) have estimated the spectral type 
of the secondary as K-M III. Since a KO III star has an ef- 
fective temperature of ~ 4800 K (e.g. Straizys & Kuriliene 
1981), we can use this as an additional constraint to limit the 
possible evolutionary histories of GRS 1915-1-105. Inspection 
of Table 1 then suggests that acceptable models for GRS 
1915+105 can have initial secondary masses as high as 6 M©. 
Indeed the model parameters for GRS 1915-1-105 in the 6 M© 
sequence are very close to the system parameters deduced 
by Greiner et al. (2001; A/bh ~ 14Mo, Md ~ 1.2 M©), for 
an assumed system inclination angle of ~ 70° , as determined 



from the orientation of the jets (Mirabel & Rodriguez 1994). 
We note that in principle one could use the 'transient' nature 
of GRS 1915-1-105 to further constrain the evolutionary past 
of the system. However, considering the unusual X-ray be- 
haviour of GRS 1915-1-105, which is very different from the 
predictions of the simple disc-instability model, this is prob- 
ably not advisable (see the discussion at the end of § 3.2). 

Our calculations have several important implications for 
GRS 1915-1-105. First, they show that the mass of the black 
hole may have increased significantly, by up to ~ 4 M© , even 
if the mass accretion rate is, on average, Eddington limited 
(also see LBW who obtained a similar estimate) . Hence the 
present mass of the black hole is not a good indicator of 
the initial black-hole mass, and any analysis of the implica- 
tions of observed black-hole masses for the evolution of the 
black-hole progenitor has to take this into account. We have 
calculated some additional evolutionary sequences starting 
with a lower black-hole mass of 7 M© (similar to the masses 
found in other black-hole binaries) and obtained acceptable 
models for GRS 1915-1-105 with black-hole masses as high 
as 11 M0. 

Second, since the black hole may have accreted a sub- 
stantial amount of matter, it may also have been spun up 
significantly and may have acquired a spin parameter as high 
as a ~ 0.8 (assuming an initially non-rotating black hole). 
This may be important for modelling the jets and the emis- 
sion from the inner parts of the accretion disc (e.g. Zhang, 
Cui & Chen 1997; Makishima et al. 2000). 

Third, the secular mass-transfer rate can be as high 
as ~ 3 X 10~^MQyr~^, implying an X-ray luminosity as 



14 Podsiadlowski, Rappaport & Han 



high as ~ 2 X 10 ergss"^. This is a factor of a few lower 
than the peak X-ray luminosity of 7 x 10"^^ ergss"^, deter- 
mined for GRS 1915-1-105 by Greiner, Morgan & Remillard 
(1996). This implies that only a moderate amount of super- 
Eddington accretion is required to explain the observed peak 
luminosity in our models to explain the observed luminosity. 

Fourth, our models predict that the surface abundance 
of the secondary should be substantially enhanced in helium 
and could show CNO abundance ratios that are close to 
the equilibrium ratios for the CNO cycles (in particular for 
the more massive secondaries). This provides a potentially 
powerful test that may help to further constrain the nature 
and the initial mass of the secondary. 

4.2 Ultraluminous X-ray sources 

Ultraluminous X-ray sources (ULX) are luminous X-ray 
sources outside the nuclei of external galaxies, typically de- 
fined to have an X-ray luminosity larger than 10"^^ ergss"^. 
They were originally discovered by Einstein (Fabbiano 1989) 
and have been found in large numbers by ROSAT and 
most recently Chandra (Colbert & Mushotzky 1999; Roberts 
& Warwick 2000; Colbert & Ptak 2002; Lira, Johnson & 
Lawrence 2002; Jeltema et al. 2002). While the physical na- 
ture of these sources has remained unclear, and indeed they 
probably form a heterogeneous class of systems (Prestwich 
et al. 2002; Roberts et al. 2002), a plausible scenario is that 
they constitute the luminous tail of the stellar-mass black- 
hole binary distribution (for recent discussions see Colbert 

6 Mushotzky 1999; King et al. 2001; King 2002; Roberts 
et al. 2002). With an inferred peak X-ray luminosity of 

7 X lO^^ergss"^ (Greiner et al. 1996), GRS 1915-^105 cer- 
tainly classifies as a typical ULX (King et al. 2001; Mirabel 
& Rodriguez 1999). The connection of ULXs with 'nor- 
mal' Galactic black-hole binaries has been strengthened by 
the determination of the black-hole mass in GRS 19154-105 
(Greiner et al. 2001), which proved that ULXs may contain 
typical stellar-mass black holes rather than a previously un- 
known class of intermediate-mass black holes of 10^ - 10"* Mq 
(as suggested by Colbert & Mushotzky 1999). 

The calculations presented in this paper strongly sup- 
port the connection of ULXs with black-hole binaries, 
in particular those with intermediate-/high-mass secon- 
daries. Figure |S| shows both the X-ray luminosity (assuming 
Eddington-limited accretion) and the potential X-ray lumi- 
nosity (assuming non-Eddington-limited accretion) of the 
binary sequences. While the Eddington-limited luminosity 
reaches a maximum of ~ 4 x 10^^ ergss"^, the potential X- 
ray luminosity can be as high as ~ 3 x lO^'^ergss"^ and 
could be even higher if the mass-transfer rate is variable 
(consistent with observations; see e.g. Lira et al. 2002 and 
the discussion at the end of § 3.2). Our calculations there- 
fore show that the mass-transfer rates in these sequences are 
high enough to provide a potential power source for ULXs. 
It requires only that the majority of ULXs have to radiate 
at a moderately super-Eddington luminosity, as is actually 
observed in GRS 1915-1-105. Even for the most luminous ob- 
served ULXs, the luminosity has to exceed the Eddington 
luminosity by a factor of ^20, which may not pose a seri- 
ous problem in radiation-pressure dominated magnetic discs 
(Begelman 2002) . This factor may be further reduced signifi- 
cantly if the radiation is beamed (King et al. 2001). Figure|K| 



shows where the most luminous systems (dark shading) are 
expected to lie in the H-R diagram of the secondary and 
in the orbital period - secondary mass plane. They indicate 
that black-hole binaries are most likely to appear as ULXs 
when they are giants or supergiants where mass transfer is 
driven by the nuclear evolution of the secondary, a phase 
that can last up to several lO'' yr (also see Fig. |S1 and Ta- 
ble 1). 

If ULXs in external galaxies are associated with 
intermediate-/high-mass black-hole binaries, they should be 
preferentially found near star-forming regions. Indeed there 
is some evidence that many ULXs are found in regions of ac- 
tive star formation, starburst galaxies and interacting galax- 
ies (Lira et al. 2002; Roberts et al. 2002; Terashima & Wil- 
son 2002; Zezas et al. 2002). On the other hand, Colbert 
& Ptak (2002) found that the number of ULXs per galaxy 
is actually higher in elliptical than in non-elliptical galaxies 
and that there is a population of ULXs in the halos of ellip- 
tical galaxies, where recent star formation is not expected 
to have occurred. This may suggest that our models are 
not directly applicable to these ULXs, unless these galax- 
ies experienced some relatively recent star formation (e.g. 
as a result of some previously unrecognized merger activ- 
ity). King (2002) has recently proposed that there are two 
classes of ULXs: relatively massive systems associated with 
young stellar populations and low-mass systems found in 
old populations. Indeed there is some evidence that a large 
fraction of ULXs in elliptical galaxies are located in globular 
clusters (see e.g. Angelini, Loewenstein & Mushotzky 2001; 
Kundu, Maccarone & Zepf 2002; White, Sarazin & Kulkarni 
2002; Jehema et al. 2002). These could be relatively long- 
period systems with giant donors (at the current epoch) that 
formed dynamically in the dense cluster cores by processes 
not considered in the present study. 

In principle we could combine the results of our evo- 
lutionary calculations with the BPS model in § 2 and pre- 
dict the properties of black-hole binaries in a typical galaxy 
(e.g. the X-ray luminosity function; see Roberts & Warwick 
2000). However, there are a large number of uncertainties in- 
volved, including the histories of stellar and globular-cluster 
formation in galaxies (especially for ellipticals), as well as 
those associated with the modelling of the evolution that 
leads to the formation of a black-hole binary and the prob- 
lem of how to relate secular mass-transfer rates to observable 
X-ray luminosities. In view of these, we cannot yet hope to 
reliably compute a population synthesis of black-hole bina- 
ries in external galaxies, and therefore little of substance is 
to be gained from such an exercise. However, putting aside 
some of these issues, e.g. the question of super-Eddington 
accretion, and taking the potential X-ray luminosity as a 
measure of the X-ray luminosity, one can obtain a qualita- 
tive idea of the relative number of systems by considering 
how much time systems spend at various mass-transfer rates. 

In Table 1 we list the X-ray lifetimes for transient phases 
and various ranges of potential X-ray luminosity in our 
binary sequences. Systems with relatively low-mass secon- 
daries are expected to be transients for most of their evolu- 
tion, while the systems with the most massive secondaries 
spend a large fraction of their X-ray lifetimes as ULXs. Con- 
sidering that our standard BPS model predicts a fairly uni- 
form distribution of secondary masses, at least above some 
characteristic minimum mass (see Figure I^J, and taking the 
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11 M© donor sequence as typical, we then estimate that 
some 10 per cent of black-hole binaries should be ULXs, 
and some 10 per cent of these have potential X-ray luminosi- 
ties above 10^° ergss"^. This would appear to be in qualita- 
tive agreement with the luminosity distribution of Roberts 
& Warwick (2000) who found that the number of sources 
in each decade of X-ray luminosity decreases according to 
L^'^ . Estimates of the absolute numbers of such binaries 
residing in a typical Milky- Way type galaxy can be made 
by taking the formation rates from Figure |3 for black-hole 
binaries with different initial donor masses and multiplying 
them by the various lifetimes given in Table 1. Such an ex- 
ercise shows that for conventional values of the A parameter 
near ~ 0.5 the predicted numbers of ULXs is substantially 
in excess of what is observed. However, for values of A closer 
to more realistic values of ~ 0.1, the computed numbers of 
ULXs is quite reasonable. 

Finally we note a piece of observational evidence that 
suggests the high luminosities inferred for some of the ULXs 
on the assumption of roughly isotropic emission are, in fact, 
correct. PakuU & Mirioni (2002) have studied very large 
(~ 300 pc) and luminous ionization nebulae surrounding 
about a dozen ULXs. They use the emission lines from these 
nebulae as an interstellar medium calorimeter to infer the 
long-term (e.g. ~ 10* yr; see e.g. Chiang & Rappaport 1996) 
mean X-ray luminosity of the underlying ULXs. From the 
structure of the ionization nebulae they can also rule out 
in some cases any extreme beaming of the radiation. The 
results of this study seem to point to good agreement be- 
tween the X-ray luminosity determined from Chandra and 
ROSAT and that found indirectly from studies of the ion- 
ization nebulae. 

4.3 Cygnus X-1 

One of the most famous black-hole binaries is Cyg X-1 with 
an orbital period of 5.6 d and a massive secondary (HDE 
226868) of spectral type 09.7 lab (Walborn 1973; Gies & 
Bolton 1986). It has a mass function of 0.252 ± 0.010 Mq 
(Gies & Bolton 1982) and used to be one of the best early 
black-hole candidates before the identification of low-mass 
black-hole transients with much larger mass functions. How- 
ever, the evolutionary state of the system has so far not been 
properly established; in particular it is not clear how much 
mass the donor star has already lost. 

In Figure|3we present the results of a binary calculation 
that may represent the evolution of Cyg X-1 and which il- 
lustrates several characteristic properties of a massive black- 
hole binary. In this model, the initial masses of the black hole 
and the secondary are 12 and 25 M©, respectively, and the 
secondary starts to fill its Roche lobe near the end of its 
main-sequences phase, when its central hydrogen mass frac- 
tion has been reduced to 0.054. The orbital period at this 
point is 6.8 d. Unlike our previous calculations we also in- 
cluded a stellar wind from the secondary of 3 x 10~^ M0 yr~^ 
(Herrero et al. 1995), taken to be constant throughout the 
evolution. 

The general evolution is reminiscent of that of an 
intermediate-mass neutron-star binary (see Fig. 7 in PRP 
and the associated discussion). After a brief turn-on phase, 
mass transfer occurs initially on the thermal timescale of 
the envelope reaching a peak mass-transfer rate of ~ 4 x 



10""^ M0 yr~^. Once the mass of the secondary has been 
reduced to a value comparable to the black hole, the sec- 
ondary reestablishes thermal equilibrium and becomes de- 
tached. Indeed because of the continuing wind mass loss 
the donor shrinks significantly below its Roche lobe dur- 
ing this phase and the system widens. The secondary starts 
to expand again after it has exhausted all of the hydrogen 
in the core and fills its Roche lobe for a second time. In 
this phase, the mass-transfer rate reaches a second peak of 
~ 4 X 10~* M0 yr~^, where mass transfer is driven by the 
evolution of the H-burning shell. The calculation was ter- 
minated at this stage, but the secondary would ultimately 
become a ~ 8M0 helium star, quite possibly becoming a 
black hole itself in the final supernova explosion. 

The most interesting feature of this calculation is that 
the system becomes detached after the initial thermal 
timescale phase because of the stellar wind from the sec- 
ondary (it acts both to widen the orbit and to shrink the 
stellar radius). During this phase, mass transfer continues 
via the stellar wind. Since the secondary is close to filling its 
Roche lobe, such a wind may be focused towards the accret- 
ing black hole, as has been inferred from the tomographic 
analysis of the mass flow in Cyg X-1 by Sowers et al. (1998). 
In this particular calculation, the secondary has a tempera- 
ture of ^ 31000 K at an orbital period of 5.6 d, in excellent 
agreement with the 09.7 spectral type of HDE 226868. It 
is worth noting that such a phase will generally not exist 
for high-mass neutron-star X-ray binaries since, because of 
the more extreme mass ratio in these systems, Roche-lobe 
overflow will generally become dynamically unstable lead- 
ing to the spiral-in of the neutron star in the envelope of 
the massive secondary. If Cyg X-1 is in the phase described 
above, the model makes the firm prediction that the sec- 
ondary should be significantly helium enriched and its sur- 
face composition should show strong evidence for CNO pro- 
cessing (in this particular model, the surface helium mass 
fraction is 0.55, i.e. roughly twice solar, at an orbital pe- 
riod of 5.6 d). Interestingly, both Herrero et al. (1995) and 
Canalizo et al. (1995) claim to have determined just such 
abundance anomalies in the secondary of Cyg X-1, which 
they argue cannot be accounted for by uncertainties in the 
atmosphere modelling. 

One caveat is that this model also predicts that the mass 
of the secondary should at the present time be comparable 
to, or lower than, the mass of the black hole. This does not 
appear to be consistent with the analyses of Gies & Bolton 
(1986), using rotational velocities and assuming synchronous 
rotation, and Gies et al. (2002), based on an Ha emission 
line analysis. Both studies suggest that the secondary is a 
factor ~ 2 - 3 times more massive than the black hole. 

However, irrespective of whether this particular model 
is applicable to Cyg X-1, the calculation in Figure |7| illus- 
trates that it is generally more likely to observe a high-mass 
black-hole X-ray binary in the relatively long-lived wind 
mass-transfer phase following the initial thermal timescale 
phase which only lasts a few 10* yr. In this example, the 
wind phase lasts a few 10^ yr, but it could last as long as a 
few lO*" yr if the secondary were initially less evolved^. This 



^ Note that, if the initial mass ratio were more extreme, the initial 
mass transfer would become dynamically unstable, as in the case 
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Figure 7. Evolutionary model for Cyg X-1. H-R diagram (top panel) and key binary parameters as a function of time since the beginning 
of mass transfer (with arbitrary offset) . Middle left: Radius of the secondary (solid curve) and Roche- lobe radius (dashed curve) . Middle 
right: Orbital period. Bottom left: Mass of the secondary (solid curve), the black hole (dashed curve) and the mass of the hydrogen- 
exhausted core (dotted curve). Bottom right: Mass-loss rate of the secondary. The secondary initially has a mass of 25 M© and has 
already exhausted most of its hydrogen in the core (it has a core hydrogen mass fraction X = 0.054). The initial mass of the black hole 
is 12 Mq and barely changes during the evolution. The calculation includes a constant stellar wind from the secondary with a mass-loss 
rate of 3 x 10~^ M© yr"'^, consistent with the present observations of the secondary in Cyg X-1. The black-hole accretion rate is always 
limited by the Eddington accretion rate. 



means that any secondary observed in a massive black-hole 
binary is likely to have already lost a significant fraction 
of its mass and that it is generally not valid to deduce the 

of massive neutron-star binaries, and the black hole would spiral 
into the massive star. In this case, there would be no subsequent 
wind mass-transfer phase. 



mass of a secondary based on its spectral type alone, as is 
frequently done in the literature. 

4.4 Low-mass black-hole binaries 

In the binary population synthesis study discussed in § 2 
and the detailed binary evolution calculations presented in 
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Figure 8. Black-hole binary formation rates as a function of the parameter fm, the fraction of the primordial primary's stellar wind 
that is ejected from the binary with the specific angular momentum of the primordial secondary. The degree of orbital separation (or 
contraction) with wind mass loss depends sensitively on this parameter. The remaining fraction (1 — /m) is assumed to be lost with the 
specific angular momentum of the primary. The value of the A-parameter has been fixed at an illustrative value of 0.1. Labels on the 
different curves refer to the maximum mass of the secondary star. Solid curves are without the inclusion of tidal interactions between 
the primary and the orbit. Dashed curves are for an assumed maximum tidal interaction, i.e. the case where the primary always rotates 
synchronously with the orbit. Note that tidal interactions typically reduce the formation rate by tending to bring systems to Roche-lobe 
overflow before the primary has lost sufficient mass to help in the process of unbinding the envelope. 



§ 3 hardly any black-hole binaries were formed with low- 
mass donors, and all the systems that were evolved in detail 
developed into wide, rather than close binaries. Thus, 9 of 
the 17 known black-hole binaries that have low-mass donors 
and Porb ^ 1 day cannot be produced within the formation 
scenario we have been evaluating, at least not without some 
modification of the input physics we have adopted. There 
are two basic reasons why such low-mass black-hole bina- 
ries do not form within the context of the model presented. 
The first has to do with the fact that the orbital energy 
in primordial binaries with low-mass secondaries is gener- 
ally insufficient to unbind the envelopes of the massive pri- 
maries (see § 2) . The second reason is that we have assumed 
that once a black-hole binary is formed with an intermediate 
mass donor, such systems do not experience magnetic brak- 
ing as a source of orbital angular momentum loss. This, in 
turn, assures that systems with initially intermediate-mass 
donors will commence Roche-lobe overfiow (RLOF) only af- 
ter undergoing a significant amount of nuclear evolution, 
subsequent to which the binary evolution necessarily leads 
to wider orbits - including periods considerably longer than 
a day. In this section we reevaluate some of our assumptions 
and relax a number of constraints to determine whether low- 
mass, compact, black-hole binaries can be formed within the 
standard scenario. 

One promising possibility, discussed in § 2, is that 



RLOF in the primordial binary might begin while the pri- 
mary has a radius between -Rhg and -Rmax, i.e. after it passes 
the HG and before it reaches its maximum radial extent. 
During this phase the star is expected to have lost a signif- 
icant fraction of its mass in a stellar wind, and this would 
greatly reduce the amount of orbital energy required to re- 
move the residual envelope in a common-envelope phase. 
However, as also discussed in § 2, the mass loss tends to 
make the orbit expand even faster than the primary can ex- 
pand, and hence RLOF is not likely to occur. In this case, 
a close binary will not be formed, if it remains bound at all. 
The calculations that led to this conclusion involved the as- 
sumptions that (1) the specific angular momentum carried 
away by the wind of the primary has the same value as that 
of the primary itself, jp, and (2) no significant synchroniz- 
ing tidal torques act between the expanding primary and 
the orbit. 

If we suppose that a fraction, /„, of the primary's stellar 
wind is deflected by the secondary and leaves the binary 
with the larger specific angular momentum of a low-mass 
secondary, js , then the orbital expansion will be diminished, 
and may even be reversed. This, in turn, could enable RLOF 
to occur when the primary has lost a significant fraction of 
its envelope. We have carried out a series of BPS calculations 
where the specific angular momentum carried away by the 
wind of the primary is set equal to 
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jw = (1 - fw)jv + fnijs- (11) 

The results for the black-hole binary formation rates vs. the 
fraction f^u are shown as solid curves in Figure m These 
were calculated for an assumed, illustrative, fixed value of 
A = 0.1. As in Figure 121 the different curves are for var- 
ious limits on the mass of the donor star in successfully 
formed black-hole binaries. Note that for = 0, the results 
agree with those in Figure |21 and show that black-hole bi- 
naries with donors < 6Mq hardly form with such a small 
value of A. However, as the value of is increased by a 
small amount, e.g. to ~ 0.05, the formation rates of systems 
with lower-mass donors grow significantly. For still larger 
values of /„,, the formation rates drop back down again. 
This can be understood as follows. For small values of /„, 
the orbital expansion can be substantially limited, and the 
primary has a chance to shed a significant part of its enve- 
lope before the common-envelope phase. For larger values of 
fw the orbit actually shrinks sufficiently rapidly that RLOF 
may commence before much envelope mass can be lost in 
a wind. The main problem with this hypothesis (invoking 
extra angular-momentum loss) is that for lower-mass secon- 
daries the fraction of the wind defiected by the secondary 
goes as the square of the mass ratio, and such fractions are 
typically much smaller than the values of fw for which this 
effect is important (however, see, Hachisu, Kato & Nomoto 
1999). 

We have carried out several other related exploratory 
tests in which the wind loss rate from the primary was taken 
to be lower than that given by the Nieuwenhuijzen & de 
Jager (1990) prescription by factors of 2 and 5. This has a 
similar effect on inhibiting the orbital widening as discussed 
above when the specific angular momentum of the wind was 
enhanced. However, for neither diminution factor of the wind 
loss rate were a significant number of low-mass black-hole 
binaries formed for any value of A. 

Finally, in this regard we tested the effects of intro- 
ducing a strong tidal coupling between the orbit and the 
rotation of the primary as it evolves to become a giant. We 
considered the extreme case where the primary remains in 
corotation with the orbit as it expands through the giant 
phase. This is equivalent to assuming that the tidal synchro- 
nization timescale is always short compared to the nuclear 
evolution timescale of the primary - which may not be un- 
reasonable for these largely convective stars (see e.g. Lecar, 
Wheeler & McKee 1976; Zahn 1977; Verbunt & Phinney 
1995). The tidal interactions have the effect of slowing the 
orbital expansion or causing orbital contraction as the pri- 
mary's moment of inertia grows, and in the case of low-mass 
secondaries the Darwin instability is even likely to set in. In 
the latter case, we simply assume that RLOF commences at 
that point. The results of our BPS calculations for strong 
tidal interactions are shown in Figure |H| as a set of dashed 
curves (same mass color coding as for the no tidal interac- 
tion case) . These represent the formation rates for black-hole 
binaries with different secondary masses - again as a func- 
tion of the /m parameter discussed above. For systems with 
higher mass secondaries, the assumed tidal interaction does 
not significantly change the formation rates. However, for 
the lower mass secondaries, it is clear that tidal interactions 
serve only to decrease the formation rates, rather than en- 
hance them. This is due to the fact that for small values of 



q = Ma/Mp, the tidal influence is quite large, leading to ei- 
ther rapid orbital shrinkage or a runaway Darwin instability. 
This forces an earlier commencement of RLOF and leaves 
much of the primary's envelope to be ejected in the common 
envelope. 

The net conclusion of these studies is that it is still quite 
problematic to form low-mass black-hole binaries directly 
via a common-envelope scenario. If, on the other hand, sys- 
tems with initially intermediate-mass donors could evolve 
to short orbital periods with low-mass donors, then initially 
low-mass donors might not be required at all. Indeed, as was 
shown in PRP, intermediate-mass X-ray binaries contain- 
ing neutron stars rather than black holes can be almost in- 
distinguishable from low-mass systems after an initial ther- 
mal timescale mass-transfer phase. Even systems with sec- 
ondaries as high as 3.5 Mq can evolve towards short orbital 
periods if the secondary is initially relatively unevolved. The 
reason is that, for intermediate-mass neutron-star binaries, 
mass transfer initially occurs from the more massive to the 
less massive component of the system which leads to a de- 
crease in the orbital period. When the secondaries develop 
convective envelopes, magnetic braking becomes the dom- 
inant mechanism to drive mass transfer, causing the orbit 
to shrink further. In the intermediate-mass black-hole bina- 
ries, the initial mass ratio is reversed and the systems tend 
to widen from the beginning. This situation would, however, 
be different if there were an additional angular-momentum 
loss mechanism operating in these systems. 

To test this possibility, we decided to relax our condi- 
tion for when orbital angular momentum losses via magnetic 
braking of the donor star are effective. In particular, we in- 
vestigated how black-hole binaries with intermediate-mass 
donor stars would evolve if magnetic braking were opera- 
ble, independent of the presence or absence of a convective 
envelope. In Figure |S| we show detailed binary evolution re- 
sults for systems with donor stars of initial mass 2, 3 and 
4M0 with the inclusion of continuous magnetic braking. 
Note that within some ~ 10^ yr of the commencement of 
mass transfer, the mass of all the donors has been reduced 
below ~ 1 Mq and that the orbital periods have decreased 
from their initial values of ~ 1/2 day through a range of 
several hours. The calculations were terminated when the 
systems reached their period minimum, which varied from 
76 min for the initially 4 M0 donor to 93 min for the ini- 
tially 2 M0 donor (the final composition in the former case 
is significantly hydrogen-depleted and has a hydrogen mass 
fraction of 0.34, while in the latter case the secondary is only 
moderately hydrogen depleted). 

We have no direct way of evaluating whether conven- 
tional ideas about magnetic braking can be overturned in 
this way or whether there might be some other angular- 
momentum loss mechanism in order to form the observed 
low-mass black hole binaries. However, it seems to us to be a 
more attractive alternative than invoking common-envelope 
scenarios that implicitly violate conservation of energy. 

4.5 Discussion and Conclusions 

In this paper we systematically explored the formation and 
evolution of black-hole binaries. Our binary calculations 
have shown that mass transfer is stable for a wide range of 
donor masses (up to about ~ 20 Mq for an initial black- 
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Figure 9. H-R diagram (top panel) and key binary parameters as a function of time since the beginning of mass transfer (with arbitrary 
offset) for binary sequences with initial masses of 2M0 (solid curves), 3 Mq (dashed curves) and 4 Mq (dot-dashed curves) where 
magnetic braking has been included even for secondaries with radiative envelopes. Middle left: Spin parameter. Middle right: Orbital 
period. Bottom left: Mass of the secondary. Bottom right: Mass-loss rate of the secondary. 



hole mass of 10 M©). After sufficient mass has been lost 
from the donor, an initially intermediate-mass donor can 
mimic a low-mass one and a high-mass donor can mimic an 
intermediate-mass one. We have shown that black holes can 
gain substantial mass from the companion even if accretion 
is Eddington-limited and that the black hole can be spun up 
significantly. This demonstrates that the present black-hole 
mass is not necessarily representative of the initial black-hole 
mass after the supernova in which it formed, an important 



fact to be taken into account when studying the implications 
of observed black-hole masses for single-star evolution. 

Our models are directly applicable to many observed 
black-hole binaries, where we particularly concentrated on 
GRS 1915-1-105 and Cyg X-1. Our calculations show that 
the initial mass of the mass donor in GRS 1915-1-105 may 
have been as high as ~ 6 Mq and the black hole may have 
accreted up to ~ 4Mq from its companion, being spun up 
in the process. The composition of the donor, in particular 
the helium abundance and the amount of CNO processing 
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it underwent, provides a potentially powerful test of the ini- 
tial mass of the donor star. V404 Cyg, a massive black-hole 
binary with an orbital period of 6.5 d and some of the best 
determined system parameters (Shahbaz et al. 1994, 1996), 
is also very well reproduced by our evolutionary sequences. 

Our calculations may also help us to understand the 
nature of ultraluminous X-ray sources (ULXs) in external 
galaxies, in particular those found in regions of active star 
formation. In our massive sequences, the systems reach po- 
tential X-ray luminosities as high as ^ 10'*^ ergss"^, compa- 
rable to the most luminous ULXs observed, when the donor 
stars become giants and the evolution is driven by the nu- 
clear evolution of the hydrogen-burning shell. 

We have also performed detailed binary population syn- 
thesis calculations which show that intermediate- and high- 
mass black-hole binaries can form at reasonable rates with 
plausible assumptions, while (initially) low-mass black-hole 
binaries apparently cannot, unless one is prepared to ac- 
cept violation of energy conservation. We have explored 
various possible solutions to reconcile this conclusion with 
the large number of low-mass, short-period black-hole bi- 
naries observed. These include: (1) The possibility that 
there are serious flaws in the stellar models of massive, 
evolved stars, or in the assumptions concerning the ejec- 
tion of the common-envelope should be considered. (2) The 
assumptions concerning the angular-momentum loss in the 
stellar wind from the primary have to be changed drasti- 
cally. (3) Low-mass, short-period systems are in fact descen- 
dants of intermediate-mass systems. This requires an addi- 
tional angular-momentum loss mechanism in systems with 
intermediate-mass radiative stars (e.g. continuous magnetic 
braking) . (4) Low- mass systems form through one of the al- 
ternative evolutionary formation channels suggested, e.g. in 
a triple scenario or out of the collapsed envelope of a massive 
star. 

Finally, we conclude that before some of these funda- 
mental issues have been resolved, the predictions of binary 
population synthesis studies of black-hole binaries have to 
be taken with considerable caution. 
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